CERN-TH. 7364/94 



TWO-LOOP QUARK SELF-ENERGY IN A NEW FORMALISM 
(I) OVERLAPPING DIVERGENCES 



Oh! 



GEORGE LEIBBRANDT^ 
Theoretical Physics Division, CERN, CH-1211 Geneva 23 



0\ • and 



JIMMY WILLIAMS^ 



O 
D 

Q 

cn ■ ABSTRACT 
(N 

psj I A new integration technique for multi-foop Feynman integrals, called the matrix 

■ method, is developed and then applied to the divergent part of the overlapping 

two- loop quark self-energy function iT, in the light-cone gauge n-A°'{x) = 0, = 
0. It is shown that the coefficient of the double-pole term is strictly local, even off 
mass-shell, while the coefficient of the single-pole term contains local as well as 
^ I nonlocal parts. On mass-shell, the single-pole part is local, of course. It is worth 

noting that the original overlapping self-energy integral reduces eventually to 10 
covariant and 38 noncovariant-gauge integrals. We were able to verify explicitly 
that the divergent parts of the 10 double covariant-gauge integrals agreed precisely 
^ ■ with those currently used to calculate radiative corrections in the Standard Model. 

5y I Our new technique is amazingly powerful, being applicable to massive and 

massless integrals alike, and capable of handling both covariant-gauge integrals 
and the more difficult noncovariant-gauge integrals. Perhaps the most important 
feature of the matrix method is the ability to execute the 4c<;-dimensional mo- 
mentum integrations in a single operation, exactly and in analytic form. The 
method works equally well for other axial- type gauges, notably the temporal 
gauge (n^ > 0) and the pure axial gauge (n^ < 0). 
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1. Introduction 



Traditional perturbation theory with its emphasis on Feynman diagrams 
continues to play a central role in quantum field theory. The success of the 
perturbative approach hinges decisively on the accurate computation of multi- 
loop Feynman integrals. Of course, there exists a vast variety of such integrals, 
but here we shall only distinguish between covariant-gauge Feynman integrals 
and noncovariant-gauge integrals. As the name implies, covariant-gauge Feynman 
integrals occur in theories quantized in a covariant gauge, such as the Fermi gauge 
or the unitary gauge. Noncovariant-gauge integrals, by contrast, arise whenever 
a noncovariant gauge is implemented, such as the powerful light-cone gauge [1] 
or the infamous Coulomb gauge [2]. While both types of integrals are known to 
require patience and a healthy respect for detail, it is also true that noncovariant 
integrals remain as popular as skunks at a garden party, as seen from the small 
fraction of higher-loop calculations in QCD and Yang-Mills theory [3-8] . 

According to the literature, multi-loop integrals, such as J (P^q J (P^k f{q, k), 
have been evaluated almost exclusively by the nested method [9-14], in which 
the four-momentum integrations are carried out sequentially (initial exponential 
parametrization of the propagators is assumed). For instance, for the double 
integral mentioned above, one first integrates over d'^^k, then over d'^^q (or con- 
versely). Dimensional regularization will be used throughout this paper, with 2u> 
denoting the dimensionality of complex space-time. We have no intention of re- 
viewing here the dominant characteristics and various idiosyncrasies of the nested 
approach, except to say that it has been used and abused with varying degrees 
of success. 

Instead, we should like to propose an alternative procedure to the nested 
method, called the matrix integration technique, in which the two momentum 
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integrals in J (P^qJ (P^k f{q,k) are integrated over 4a;-dimensional space in a 
single operation. Perhaps the most appeahng feature of this method is the abihty 
to perform the momentum integration exactly and in closed form, thereby guar- 
anteeing from the outset a certain amount of what might be called "calculational 
streamlining" . In the course of our investigation we shall encounter several ex- 
amples of this "streamlining" . It turns out that the matrix technique works for 
covariant and noncovariant gauges alike, and regardless whether the integrals are 
massive or massless. 

The purpose of our integration program may now be stated as follows: 

1. To develop for two-loop integrals an alternative approach to the nested 
method, called the matrix integration technique. 

2. To apply this technique to the complete two-loop quark self-energy function 
in the light-cone gauge n-A°-{x) = 0, = 0. 

3. To derive, if possible, the necessary counterterms and use these to renormal- 
ize the theory to two-loop order. 

4. To extend the matrix technique to other axial gauges, notably to the tem- 
poral gauge defined by n-A°'{x) — 0, > 0, and to the pure axial gauge 
n-A'^{x) = 0, < 0. 

As indicated above, the testing ground for our matrix technique will be the quark 
self-energy function to two loops, depicted graphically in Figures 1 and 2. For 
pedagogical reasons, we have decided to report our results in two separate papers. 
In paper I, we shall discuss various mathematical tools and then apply these to 
compute the divergent portion of the overlapping quark self-energy (Fig. 1). The 
remaining two- loop graphs, including the non-trivial rainbow diagram (Fig. 2(a)), 
will be treated in paper II, where we shall also study the counterterms required 
for renormalization. 
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The plan of paper I is as follows. In Section 2 we define the structure of 
the overlapping quark self-energy function and review some standard formulas 
needed for one-loop calculations. In Section 3, we summarize the main features 
of the matrix method, including reduction of the integrand to Gaussian form and 
the evaluation of 4a;-dimensional Gaussian integrals. The parameter integrations 
are carried out in Sections 4 and 5. But in order to prepare the reader for 
certain technical subtleties, we shall first make a short detour to introduce some 
terminology. 

For each factor in the denominator of the integrand of J (P^q J (P^k f{q, k), 
there will be a Schwinger parameter, aj say (j = 1,2, . . .), with an infinite domain 
of integration; i.e., aj e [0,oo]. For the overlapping diagram discussed in this 
article, we transform the set of parameters {aj} to a more "user-friendly" set S, 
containing two types of parameters: one type / parameter (A), with an infinite 
domain, and up to six type //parameters {X, P,G,b, h, a), with finite domains. 
We shall demonstrate in Section 4 that integration over the lone type I parameter 
leads to a simple pole, and in Section 5 that one more simple pole (and hence 
a double pole overall) emerges from integration over one of the finite type II 
parameters. The paper concludes in Section 6 with a short discussion. 

2. Basic Tools 
(a) Notation and review of the light-cone prescription 

In Yang- Mills theory, the light-cone gauge is characterized by 

n'^A''^{x) = 0, n2 = 0, // = 0,1,2,3, (2.1) 

where — (no, n) defines a fixed axis in four-space [1]. We use a (-I-, — , — , — ) 
metric and employ dimensional regularization in a space-time of 2u> dimensions. 
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The relevant SU(3)c color Lagrangian density reads [15] 



£ = lim ^2^Q^ ^2.2) 

A — ^'O 

Ant = -lF-,F->^^ + J2 #a(7M^^^ + , (2-3) 

k 

where ipa and represent fermion and gluon fields, respectively, 7^ are 4x4 
Dirac gamma-matrices, a = 1, 2, ... ,8 is the group index, a, P are color indices, 
A is the gauge parameter, k — u, d, . . . is the quark flavor index, and are 
quark rest masses. Moreover, 

F;, = d^Al - d^A; + gr^'^A^'^Al (2.4a) 

is the field strength {g is the QCD coupling constant), and 

Da(3,^ = 5^pd^-igT^pAl (2.4b) 

denotes the covariant derivative {d^ = d/dx^). T^^ are the generators of the 
gauge group SU(3) which obey the commutation relations 

[T",T^] = i/^'^^T^ (2.4c) 

jahc i^ging antisymmetric structure constants. In the light-cone gauge, the La- 
grangian density (2.2) leads to the gluon propagator 

«;t(.) = ^(..-^^^), e>0. (.5) 

The spurious poles of the factor {q-n)~^ will be handled by the n* -prescription 
[16-17]: 

— = lim , e > 0, (2.6) 

q-n e-»o q-nq-n*+ie 

with = (no,n) and n* = (no, —n). 
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(b) The overlapping quark self-energy 



Application of the Feynman rules for QCD, with the quark- quark-gluon ver- 
tex factor igT^p^'^, the quark propagator — m), and the gluon propagator 
(2.5), yields the following expression for the two-loop overlapping quark self- 
energy function in the light-cone gauge (Fig. 1): 

lo = ig^'Tl.^T^^Tl.cnj,!, (2.7) 
J f f d'^Q.d'^k 7^ {J> — ^ + m) 7"^ {j> — }^ — ^ + m) '-y^ — Ijt + m) Y 



M J M (27r)^ [{p — qY — rn^] [{p — k — g)^ — m?] [{p — k)'^ — rn?] fc^ 

^ i^-^ ) ^^ir- J ' ^'-'^ 

where m is the quark mass, p the quark 4-momentum, denotes integration 
over Minkowski space, and where the five ie terms have been omitted. Evaluation 
of the integral / by means of the matrix method is the main objective of this 
article. 

To begin with, we observe that / in Eq. (2.8) diverges at the following three 
types of limits: 

(i) \q\ and/or \k\ — > 00 ; 

(ii) q and/or k — > zeros of the quadratic factors in the denominator of Eq. 
(2.8); 

(iii) q-n — > and/or k-n — > 0. 

The first two types of divergence are handled by dimensional regularization and 
Wick rotation. The third type of divergence (iii) is symptomatic of axial-type 
gauges, and requires application of the n*-prescription (2.6) to the spurious poles 
of {q-n)~^ and {k-n)~^. Since the n* -prescription permits a Wick rotation 
(zg4 = qo, ik/^ = ko, ipA = po), we may Wick- rotate the entire integrand from 
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Minkowski space to Euclidean space. 

Next we apply exponential parametrization to every quadratic factor Fj 
(j = 1,2,...) in the denominator of the rotated integrand: 



thereby replacing the product of quadratic factors in the denominator by a sum 
of quadratic terms in the exponent. Completing the square in this exponent, 
we can then integrate over the momentum variables by means of the generalized 
Gaussian formula: 



(c) Two-loop integrations: the traditional approach 

Although the methods described in Section 2(b) are applicable to any multi- 
loop integral, the ever-increasing technical difficulties have drastically restricted 
the number of explicit computations. Here is a brief look at some typical prob- 
lems, encountered already in double integrals such as J dP^q J d'^^k f{q,k). 

Two-loop integrals have traditionally been computed by the so-called nested 
method. As alluded to in the Introduction , the basic idea behind the nested 
method is to integrate over the loop momenta one loop at a time, i.e., to integrate 
over k as completely as possible, before attempting integration over q. In many 
two-loop cases the A;-integral is so complicated, however, that we cannot fully 
compute it, unless we express the result in the form of a series, such as the Laurent 
series {u> — 2y Tj{n,p, q; m). Each of the functions Tj must then be multiplied 
by the remaining qf-dependent factors in the original integrand, and the resulting 




F,>0, 



(2.9) 




A > 0. 



(2.10) 
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expressions integrated over 2a;-dimensional q-spa.ce. One hopes, of course, that 
only the first few (i.e., small j) terms of the Laurent series will contribute to the 
pole parts of the final result. Unfortunately, no such luck prevails in general, as 
may be seen from the following typical, albeit simplified, example: 

^"gi-- ^2^-4 ^ (Ing)^ dq^ . (2.11) 

Integration of the j-th term on the right-hand side by parts j times yields Ylij 2-' , 
whereas direct evaluation of the left-hand side yields {2 — u))~^. The inconsistency 
stems from the fact that the exponential series for does not converge 

uniformly as \q\ — > oo. 

In order to improve the behaviour for large values of we could try to 
express the result of the A;-integration as a series in powers of l/q^. One may 
always obtain such a series by applying the binomial formula to the integrand 
before integrating over the last one or two Schwinger parameters associated with 
the A;-integration. On the surface, this strategy looks promising, but on closer 
inspection we sometimes find that the resulting series fails to converge for certain 
combinations of values of q and the unintegrated Schwinger parameters. Some 
of the UV divergences of the integral /, for instance, occur precisely in regions 
where the binomial series diverges. 

These simple examples serve as a potent reminder that (1) the divergent 
terms for overlapping loops do not arise merely from the divergent and finite 
terms (j < 0) of the individual loops, and that (2) care must be exercised in 
choosing series which converge uniformly over the entire region of integration. 
We shall see in Sections 4 and 5 that point (2) is also crucial in the context of the 
matrix method, where similar convergence issues arise, albeit in vastly simplified 
form. 

7 



3. The Matrix Method 
(a) Basics 



In this section we propose an alternative procedure to the nested method, 
called the matrix integration technique, or matrix method for short. In the matrix 
method we regard q-space and fc-space as subspaces of a single 4a;-dimensional 
momentum space and proceed as in the case of one-loop integrals. In other words, 
we Wick-rotate qo and ko, use the n*-prescription (2.6) for {q-n)~^ and {k-n)~^, 
and apply the exponential parametrization (2.9) to all denominator factors. The 
two-loop integral / in Eq. (2.8) for the overlapping diagram in Figure 1 then 
assumes the form: 

p f /"OO /"OO 

d^^q d^^k f{q,k) = dai da2 . . . J[P{q,k)] = , (3.1) 

Je Je Jo Jo 

where: 

i^(— 1)^ 7^^ {i) — ^ — m) Y {i> — ^ — i — ^ Y {j> — ^ — m) 



f{q,k) = 



(27r)4'^ q^ [{p - qY + w?] [{p - k - qf + w?] [{p - kf + w?] k'^ 



X- —- ; (m this case); (3.2) 

n-q n-k 



3) 



J[P(g,fc)] = / d^^q I d^'^kP{q,k) exp(-z^Mz + 2B^z - C) ; (3. 
J E Je 

P{q,k) = the numerator of /, multiplied by —q-n*/n^ and/or 
— k-n*/nQ from the n* -prescription (as applicable); 

z = {k^, q4, ks, qfs, ...)"*" = a 4a;-component column-vector ("*" = transpose); 
M is a 4a; X 4a; matrix; B is a 4a;-component column- vector; C is a scalar; 
and denotes integration over Euclidean space. Since My is the coefficient of 
ZiZj in the exponent of the integrand in Eq. (3.3), we may always symmetrize M 



C'^Zi" refers to the i^^ component of z). 

In the case of the integral /, with / given by Eq. (3.2), let us take ai, ... ,0:7 
to be the parameters corresponding to the factors n-q, , [(p — g)^ + m^], [{p — 
A;-q')^+m^], [(p-A;)^+m^], A;^, and n-k, respectively. We then find that M, B, 
and C take particularly simple forms when we change the variables of integration 
from {cKi, . . . , CKy} to the more "user- friendly" set S = {A; A, /?, G, b, h, a}, where 



(3.4) 



A= cti -I- Q!2 + 0:3 + 2a4 + + aQ + ar , 
G= a4,/A, b= (q!4 + 05)/^ , 
P= (as + a4)/A, h= (0:4 + 05 + ae)/^ , 
A= {a2 + as + a4)/A, a= (0:4 + as + cie + ci7)M • 

As explained in the Introduction , A is a type I parameter with an infinite domain, 
while A, P, G, b, h, a are type II parameters with finite domains. In terms of these 
parameters. 



M 





-a G 










- bp4 ' 




G 1-a 










(3p4 




a G 










bps 




G 1-a 








, B = A 


PP3 


A 












h 


G 








bp2 




G 


A 








PP2 








h 


G 




bpi 








G 


A. 




iPpil 



(3.5) 



C^A{b + P-G){p^ + m^), 

where ni = n2 — for simplicity. Note that the lower half of M has 2ij — 2 
pairs of rows, and the right-hand half 2a; — 2 pairs of columns, so that 



Vdet(M) = yl^^Dii 
where D» = a{l — a) — G"^ , 



D± = Xh-G'^. 



(3.6) 
(3.7) 
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(b) Momentum Integration 



Our first major task is to evaluate the momentum-space integral J[P] in Eq. 
(3.3), P being a polynomial in the components of q and k. Differentiating Eq. 
(3.3) partially with respect to Bi , we obtain 

^ = 2J[.,P]. (3.8) 

Once J[l] is known, we can derive J[P] for any polynomial P by repeated appli- 
cation of formula (3.8). 

To evaluate we first diagonalize the quadratic form in the exponent of 
the integrand in Eq. (3.3). Since M is symmetric, there exists a matrix R such 
that R"^R = 1 , and RMR"^ is diagonal. Letting z = YCy + M-^B in Eq. 
(3.3), and choosing P = 1, we obtain 



J[l] = exp(B^M-^B - C) I exp (-y^RMR^y) d^'^y 

J E 



(3.9) 



With RMR^ diagonal, the above integral is just a product of one-dimensional 
Gaussian integrals, so that 

_ .^-exp(B^M-^B-C) 
^ ^ " Vd^ ' 

det(R) and det(R^) having cancelled each other. Applying formula (3.8) repeat- 
edly to Eq. (3.10), we get 

J[z,] = miJ[l], (3.11a) 
J[ziZj\ = [rmmj + J[l] , (3.11b) 

J[ziZjZk] = [minijmk + mi{^^^)jk + TOj(^V")»fe + "^fe(^V")y] "^[^i ' (3-llc) 



where m = M-^B, and dmi/dBj = (M-i)^ . 
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When M and B are given by Eqs. (3.5), D\\ and by Eqs. (3.7), and p\\ 

and p±^ by 

„ , n*-pn + n-pn* ^ , 

Pll ^ (0, 0, p3, P4) = ' (3-12a) 

P± = (Pi, P2, 0, 0) =p-p||, (3.12b) 

we find m = (r4, S4, rs, S3, .. .)"■", (3.13) 

/(l-a)6-G'/3\ {\h-G(3\ , , , 

with r= rii+r^ = (^^ ^ j P|| + )^^' ^ 



. = .II + .± = ( ) + ( ^^^^ ) ■ (3-i4b) 



(For the sake of clarity, we have dropped the Lorentz indices on py, p±, n, n*, 
T^i T^Wi i^±i s, S||, and s±.) Substituting from Eqs. (3.5) and (3.13) into Eqs. (3.10) 
and (3.11), and exploiting the linearity of J[P], we obtain the useful momentum- 
space integrals given in Appendix A. 

The off-diagonal entries in M originate from the q-k term in the exponent 
of the integrand in Eq. (3.3). In the integral /, this term comes, of course, from 
the parametrization of the denominator factor [{p — k — qY + w^], which in turn 
comes from the internal line shared by the two loops in Figure 1. Note that in 
the absence of such a shared line, the integral I would factor into two separate 
one-loop integrals. By diagonalizing M, it might appear that we have effectively 
disentangled the overlapping loops, but that is not the case. All we have done is 
shift the problem to the parameter integrations. To see this, we observe that the 
integrands of these integrations still contain the off-diagonal elements of M by 



virtue of the factors M ^ and y/ det(M) in Eqs. (3.10) and (3.11). Parameter 
integrations will be discussed in Sections 4 and 5. 
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(c) Reduction of the integrand 



We observe in Eqs. (3.11) that the complexity of J increases dramaticaUy 
with the degree of its argument P, P being the numerator of the integrand / in 
Eq. (3.1) multiphed by —q-n*/n'Q and/or —k-n*/nQ from the n*-prescription if 
apphcable. It is desirable, therefore, to reduce the degree of P as much as possible 
before integration. For /, the initial degree of P is 7, but we can reduce this 
number to 3 by executing the following trivial operations before Wick-rotation: 

1. We completely expand the numerator of the integrand in Eq. (2.8) into a 
sum of products. We then drop all terms which are shown by power counting 
to be UV-convergent, since we are only interested in the divergent terms. We 
recall that the n* -prescription is consistent with power counting [1] , and that 
Weinberg's Theorem asserts that a two-loop integral is UV-convergent if the 
net power of the integrand in each of q, k, and both q and k together, is 
less than zero [18]. For /, the result is that terms containing more than one 
factor of m and/or ^ may be dropped. 

2. In the remaining terms, we rearrange the non-commuting Dirac matrix 
factors with the help of the relation [y^^,^'^]+ = 2gi^^ , so that most terms 
with more than one matrix factor cancel with one another. Some terms 
involving j^^i^ remain, but in these terms we may replace by + = 
q-k because of the symmetry between q and k in Eq. (2.8). Also note that 
7'^7^ = 2a;, and that = in the light-cone gauge. 

3. In each of the remaining terms, we cancel as many factors as possible with 
factors in the denominator of the integrand. (In some cases, identities such as 
2q-k — —w? — [{p — qY — m^] -\-[{p — k — g)^ — w?] — [{p — k)"^ — m^] , 
2p-q = p^ — — [{p — qY — rn?] + q^ , etc. are helpful.) The resulting 
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integrands will have fewer denominator factors than were present originally 
in Eq. (2.8), but Eqs. (3.5) to (3.14) are still valid, provided we set the a 
parameters which correspond to absent factors to zero in Eqs. (3.4), and 
omit the integrations over these parameters from Eq. (3.1). (We see from 
Eq. (2.9) that this procedure is equivalent to inserting factors of exp(O) in the 
integrand in Eq. (3.3). The advantage of this strategy is obvious: it allows us 
to use one common expression for J for all integrals with the same P, even 
though the integrands of all these integrals will have different denominators.) 

4. Finally, we reduce the number of independent terms by interchanging q and 
k in selected terms. 

The above four steps were carried out by means of a computer program 
designed by one of the authors (J.W.). Application of these steps, followed by 
Wick-rotation, transforms the integral / into the sum of 53 integrals, represented 
by the 53 non- blank entries in Tables 1 to 6. Each of the 53 integrals is equal to 
an entry in one of the tables times the Euclidean-vector expression at the left- 
hand side of its row of the table, divided by (27r)^'^ times the denominator at 
the bottom of the table, and integrated over 4a;-dimensional Qf/c-space. 

The symbols n, q, Q, a, and k at the bottoms of the tables represent 
denominator factors in accordance with the following scheme: 

n q Q K k n 

n-q q^ [(p - 9)^ + m^] \{p - k - q)^ -\- 'm^\ [(p - A;)^ + m^] k^ n-k 

(3.15) 

Thus, an "n" on the left side of the denominator represents n-q, while an "n" 
on the right represents n-k. For example, the entry 4 — 4e in the fifth row and 
second-last column of Table 1 corresponds to the Euclidean-space integral 
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(4-4e) n-p f f i (P^qcP^k ^^^^^ 



where e = 2 — uj ; e occurs in the integrand only because 7^*7^^ = 2a; (see step 
2 above). 

The zeros and ones at the right-hand sides of some table entries are explained 
in the Key following Table 1. The two entries marked 00 in Table 2 may be 
ignored, since it can be demonstrated by the nested method (without the use 
of series) that the corresponding integrals cancel each other exactly. The three 
stars (*) in Tables 4 and 6 indicate integrals whose contributions to the divergent 
parts of / vanish due to symmetric integration. Hence only 48 integrals require 
evaluation. All necessary momentum-space integrals J[P] are given in Appendix 
A, leaving only the integrations over the parameters still to be done. 



Table 1 : Terms of / with 5 denominator factors including q-n . 



1^ (p^+m^) 






2 + 4e 


-2 + 2e 


4-4e ° 


mn-k 


8-8e 


4 






8-8e [] 


-8 + 8e 


n-k 


8-8e 


-4 






8-8e [] 


-8 + 8e 


n-p ^ 


-4 + 4e 


20 -8e 






-4 + 4e [] 


4-4e ° 


n-pi 


4-4e 1 


-4 1 


4 + 8e g 


4-46 


-4 + 4e 


tji p-q 


1 


-4 1 












01 


fji k-q 


1 


1 










-4 + 46 


n-k ^ 


4e 1 


1 


4-46 


-8 + 8e [j 


4-46 


n-k Ijt 


-4 + 4e 


-4 + 4e 




00 




00 




00 




nq ^Kk 


n Q^Kk 


nqQAK 


nqQA k 


nqQ Kk 
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Key to a typical entry : 

j- value of a at which / J\ diverges. (See Section 4.) 

I 4 _ 4g I ^ 

I values of a at which J Jq diverges, oo indicates that 

I the degree is -1 at a = 0. (See Section 4(b).) 

^ — coefficient of the term ( e = 2 — a; ) . 



Table 2 : Terms of / with 4 denominator factors including q-n . 



i 


2e 1 


-2e 01 


-2-2e 


-4-2e 


-2+2e 00 


2-2e 00 




n AKk 


n Q Kk 


n Qa k 


nq aK 


nqQA 


nqQ K 



Table 3 : Covariant-gauge terms of /. 



m 


2-2e+2e2 i 


1 


8 -8e2 


2-2e+2e2 







-2+2e-2e2 


01 




4-2e+2e2 i 


1 


8-8e-8e2 o 


4-2e+2e2 







-4+2e-2e2 


01 


i 


4-4e 1 


-8 + 8e 1 




8+4e+4e2 o 


-8 +4e2 







4-4e 


01 




qAKk 


QAKk 




qQaK 


qQA k 






qQ Kk 




Table 4: Terms of / with 5 denominator factors including 


q-n 


and k-n. 






4 


-4 


8 * 


-4 


* 










nqQ K n 


nqQA n 


n QaK n 


nq aK n 









Table 5 Table 6 



n-k Ij^ {p^ + w?) 


4-4e ° 


i/i n-p (q-k)^ 


4 * 




nqQAKk 




nqQAKkn 
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4. Finding the Divergences 
(a) Integration over infinite-parameter space 



In order to complete the parameter change from {aj} to {A; A, (3, G, b, h, a}, 
we have to apply the transformation (3.4) to the integrations on the right-hand 
side of Eq. (3.1) for each of our 48 integrals. Integrals containing all seven "fac- 
tors" n, q, Q, A, K, k, n in their denominators then transform as follows: 



nOO nOO 

IeU] = da^... darJ[P] 
Jo Jo 



pl — G pl—Q, M i>X i>h poo 

dG da dX dh / d(3 db A^dAJ[P]. (4.1) 
Jo Jg Jg Jg Jg Jg Jo 



But let's suppose now that a certain denominator factor is missing from the 
simplified integrand / in Eq. (3.1). In that mentioned in Section 3(c), 

we must set the corresponding a parameter to zero in Eqs. (3.4), and omit 
integration over this parameter in Eq. (3.1). To see the effect of these changes 
on the transformed integral (4.1), consider the example (3.16), in which / = 
^/(nqQ/^k). Since the fifth and seventh denominator factors are absent in this 
case, we put = = to get 



lE[f] = [ d^^q [ d^^k 
Je J e 



nqQA k 



/■oo pOO pOO pOO fOO 

= / dai / da2 / das / da4 / daQ J 
Jo Jo Jo Jo Jo 



-n*-q ^ 



Hr 



CUs = Q!7 = 



/■i f-l-G pi-a rX 

dG da dX dp 
Jo Jg Jg Jg 



Ja{X,(3,G, b,h,a) 



b^G 
h ^ a 



(4.2) 



with Ja 



Jo 



AUAJ 



—n*-q 



POO 44 

_ / dA — rJ[n*-qq^,]. (4.3) 

^0 "-0 
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In the next few pages, we shall concentrate on the detailed evaluation of the 
integral (4.2). 

In view of Eqs. (A. 4) and (A. 8), we find that J a in Eq. (4.3) becomes 
il Jo D\\DY 



Integration over A is straightforward. Note that Eq. (4.4) shows the entire A- 
dependence of the integrand, since r, s, Dy, Dj^, and H have been defined so 
as to be independent of A. Using Eqs. (3.14), together with 



poo 

j t^e-*dt = r(a; + l), Re(x + 1)>0, (4.5) 
Jo 



we find that 



Ja 



7^4-26 



h ^ a 



[r(2e) Jo + r(l + 2e) Ji + r(2 + 2e) J2 ] , (4.6) 



where 



and J2 = 0. Similar results may be obtained for the other integrals in Tables 1 
to 6; only the expressions for Jq, Ji, and J2 vary. 

As explained in the Introduction , the first pole term r(2e) Jo in Eq. (4.6) 
arises from integration over the type I parameter A. This divergence may be 
traced back to the fact that the original integral / is UV-divergent with respect 
to q and k taken together. Additional UV subdivergences will emerge when the 
Jo and Ji terms are integrated over the finite parameters A, /3, G, b, h, and a 
(called type II parameters in the Introduction). These subdivergences arise from 
the two loops taken separately. 
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In view of the appearance of , D_i, and H in Eq. (4.7), exact integration of 
Jo and Ji over the finite parameters would pose a major challenge. Fortunately, 
however, we are only interested in contributions to / which diverge as a; — > 2; 
i.e., in the terms with negative indices in the Laurent series 

/(e) = /_t(p,n;m)e-* + /i_t(p, n; m)e^-^ + ... , (4.8) 

where t is an integer which, in our case, happens to be 2. The divergences occur 
only at certain boundaries of the region of integration in parameter space. We can 
extract the divergent terms exactly by integrating the first one or two terms of 
suitably chosen series which converge to J[P] near these boundaries. In choosing 
these series, we must be careful to avoid the types of problems which plague the 
nested method. 

In the case of the Jq term, one boundary at which the parameter integrations 
diverge as a; ^ 2 is just A = 0; this fact explains why integration over A 
generates the divergent factor r(2e) in this term. Of course, the presence of this 
factor implies that the coefficient /_i in Eq. (4.8) will depend on the finite part 
of the integral of Jq over the remaining parameters. It would seem, therefore, 
that the Jq terms will have to be integrated over the entire region of integration 
$ of these parameters. 

Fortunately, things are not quite as bad as they seem, since the most com- 
plicated part of the Jq term is the factor H~'^^. If we use the expansion 

H-^^ = e-2^^n^ = 1 - 2e\nH + 2e^{\nHf - ... , (4.9) 

it would appear that we only need to keep the first term of the series when inte- 
grating over the whole region since the common factor of 2e in the other terms 
will cancel the divergence arising from the factor r(2e) in Eq. (4.6). However, it 
was just this kind of reasoning that led to the catastrophe in the nested method 
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exemplified by Eq. (2.11). To avoid a similar catastrophe here, we must verify 
that Jo is convergent at all zeros and poles of H before we may use the series 
(4.9). It follows from definition (A. 9) that H has no poles in but H does go to 
zero when 6 + /? — G— >0. In order to determine whether the integral of any Jq 
term diverges at any of these zeros as co ^ 2, and in order to enable us to find 
all divergent contributions to the final result /, we must systematically determine 
the location in finite-parameter space of every subdivergence of every Jq, Ji, and 
J2 integral. This task will be tackled in the next subsection. 



(b) Subdivergences 

Divergences in the finite-parameter integrations, as a; — > 2, could potentially 
occur at the zeros of any of the three factors -Dy , D±, and H which appear in the 
denominators of the Jq, Ji, and J2 terms. These zeros occur on certain bound- 
aries of the integration region such as a = 0, X = h = G, etc. Specifically, 

D|| — > linearly if a — > 0, a — > 1, or a,G ^ ^ , 

quadratically if both X, h ^ , } (4.10a) 

or else linearly if A — > 0, h ^ 0, or \,h ^ G , 



(cf. Eqs. (3.7) and (4.1)). To locate the zeros of H, we observe that the exponent 
of the integrand in Eq. (3.3) is less than or equal to zero by construction, and 
equal to -AH when y = (cf. Eqs. (A.9) and (3.9)). Hence H >{) even if the 
mass m vanishes, while for m ^ 0, 

H ^0 linearly, if and only if both b, (3 ^ . (4.10b) 



Depending on the specific form of the integrand, a given finite-parameter 
integral could conceivably diverge at the intersection of some particular subset 
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of the above boundaries, but not at any other points. Accordingly, we must 
consider all possible intersections of these boundaries for every term of every 
finite-parameter integrand. The set of possibilities which needs to be examined is 
severely constrained by two chains of inequalities, which follow from Eqs. (3.4): 

a>h>b>G>0 and l-a>A>/3>G'>0. (4.11) 

If any particular parameter in either of these chains goes to zero, then all param- 
eters to its right in the chain must also approach zero. 

In order to determine whether one of our integrals is divergent at a particular 
boundary of we observe that the integral of a rational function of the param- 
eters diverges at some point in finite-parameter space only if the degree of the 
integrand (including the measure dG da . . .) is less than or equal to zero at this 
point. The degree of the integrand may be calculated according to the following 
rules: 

1) A linear function of a parameter has degree 1 if it goes to zero at the point 
in question; otherwise the function has degree 0. 

2) The degree of a product (quotient) is the sum (difference) of the degrees of 
its parts. 

3) The degree of a sum or difference, in the numerator, is the minimum of the 
degrees of its parts. For example, if / ^ linearly and g quadratically, 
then f-\-g^{) linearly. (If the parts have equal degree, this rule may predict 
divergences which eventually cancel each other.) 

4) The degrees of the denominator factors Dy, and H are given by Eqs. 
(4.10) for all boundaries except those at which the degrees are zero. 

5) The degree of the measure is the dimensionality of finite-parameter space 
minus the dimensionality of the boundary in question (not counting dimen- 
sions associated with absent parameters, such as h and h in Example (4.2)). 
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For instance, at the boundary a = , inequalities (4.11) imply that G must 
also be zero; hence the degree of the measure dG da dX df5 at this boundary 
is 4 — 2 = 2 . To see why, consider a transformation of a and G to polar 
co-ordinates. 

As an example, consider the three terms comprising Jq and Ji in Eqs. (4.7). 
If we multiply these terms by dG da dX d(5 to complete the integrands, we find 
that the degrees of these integrands are as shown in Tables 7 and 8. Each number 
in these tables represents the degree of the integrand at the boundary at which 
the parameters, shown above it and to the left, go to zero. From these tables we 
see that all three terms in Eqs. (4.7) have subdivergences at the boundary where 
a — > but not /? — > or A — > 0. 

Table 7 : Degree of Jq Table 8 : Degrees of Ji terms 

at various boundaries. at various boundaries. 





G 


/3 


A 1 -a 




G 


/3 


A 1 -a 




G 


/9 


A 1 -a 


G 


1 


2 


2 1 


G 


1 


3 


3 1 


G 


1 


3 


2 1 


a 





1 


1 


a 





2 


2 


a 





2 


1 



When tables, similar to Tables 7 and 8, are constructed for every term of 
every finite-parameter integral which arises in the evaluation of /, it is found that 
subdivergences occur only at the following two boundaries of 

1) a — > and not P ^ (corresponding to divergent k integration). 

2) a — > 1 and not 6 — > (corresponding to divergent q integration). 

Each integral which diverges at one or both of these boundaries is marked in 
Tables 1 to 6 with ° and/or ^, as explained in the Key. Note that neither of 
these types of divergence occurs at a boundary where b and /? both go to zero. 
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In particular, we see from Eq. (4.10b) that there are no subdivergences at zeros 
of H. Therefore, it is safe to use the series (4.9) in aU cases, as we shaU do from 
now on. 

To calculate the divergent contributions to / which arise at either of the 
boundaries a = or a = 1 , we expand each integrand J, (i = 0, 1, 2) as a series 
in powers of the parameters which go to zero at that boundary, and keep only 
those terms whose integrals diverge. (We justify this procedure by using the five 
rules given above to show that the integral of the difference between Jj and the 
terms we keep is convergent at all boundaries in every case.) In all but two of 
the cases, the degree of the original integrand is zero at the subdivergence, so 
that only the leading term of the series is needed. For each of the other two cases 
(marked oo in Table 2), the degree of the Jq integrand is —1, and two terms 
would be needed. Fortunately, we are saved from having to integrate these terms, 
since these two exceptional integrals cancel each other exactly (cf. Section 3(c)). 

5. Finite-Parameter Integration 
(a) Integrations at subdivergences 

Employing the two series expansions discussed in Section 4, we can now 
complete the finite-parameter integrations to the extent necessary to obtain ex- 
act expressions for the coefficients /_i and /_2 in Eq. (4.8). Proceeding with 
Example (4.2), we first substitute the series for H~'^^ in Eq. (4.9) into Eqs. (4.7). 
We then use Eq. (4.6) and the identity T{z + 1) = zT{z) to obtain 



TT 



.4-2e 



Ja 



b 

^ h 




G ~ 



n. 



2 




[r(2e)yo + r(l + 2e)yi + r(2 + 2e)y2] +0(e). 



a 



(5.1) 



22 



n = 




) 



1^* a InH 



e ' 



(5.2b) 



and Y2 = 0. Similar results may be obtained for the other integrals in Tables 1 
to 6; only the expressions for Yq^ Yi, and I2 vary. It remains to integrate these 
expressions over the finite parameters. Since Yi and Y2 are not multiplied by 
divergent F-functions, they contribute to /_i only through the subdivergences 
in their finite-parameter integrals. The single Y2 integral has no subdivergences, 
and may be ignored. The 0{e) term in Eq. (5.1) may also be ignored, since there 
are no subdivergences at the zeros of as explained in Section 4(b). 

Returning to our example in Eqs. (5.2), we now compute the divergent part 
of the integral of over the finite parameters (3, A, a and G (cf. Eq. (4.2)). We 
recall from Section 4(b) that all subdivergences of Eq. (4.2) occur at the boundary 
a = G = 0. Therefore, we expand Yi as a series in the parameters a and G, and 
integrate only the divergent leading term. To derive this term, we simply replace 
H by its leading term Hq = Ha^o , and drop and (1 — a) from D\^, D±, and 
from the numerator in Eq. (5.2b). In this way we obtain 



Applying Eqs. (A. 9), (3.14), and (3.7), and recalling that a = h in our example, 
and that < 6 < a (cf. inequality (4.11)), we find 



Yi = E + 0{a 




(5.3) 



where 




(5.4) 




(5.5) 
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Integration of Eq. (5.3) over /5, A, a, and G in accordance with Eq. (4.2) yields 



Yi = I dG I da dX dp a^~^ E + finite terms , 



* JO JG JG JG 

1-G 



= pdG I 

Jo JG 



a^ ^ da (Lq — Li) + finite terms , (5.6) 

G 

where Lq = [ dX [ dp E , (5.7a) 
Jo Jo 

rl rX pG rl 

Li = dX dp E + dp dXE, (5.7b) 

Ji-a JG Jo Jp 



and E is defined by Eq. (5.4). Integral (5.7a) is convergent even for a; = 2 (e = 0), 
and may be evaluated in the e = case with the help of formulas (B.9), (B.8), and 
(B.5) from Appendix B. The two integration regions in Eq. (5.7b) are entirely 
inside the region in Eq. (5.7a), and shrink in proportion to a. Consequently, 
Li^Oas a^O; hence the integral of Li over a and G may be absorbed into 
the finite terms in Eq. (5.6). Lq may then be factored out of the integral, since it 
is independent of a and G. The a and G integrations are easily completed: they 
produce the divergent factor e~^, plus a finite term, so that Eq. (5.6) becomes 

f Yi = - f dX f dp Ee=o + finite terms . (5.8) 
i$ £ Jo Jo 



The Yi terms arising from the other integrals in Tables 1 to 6 may be in- 
tegrated in similar fashion. After convergent terms have been discarded, each 
integral factors into a trivial divergent integral, times a finite integral which can 
be evaluated at e = with the help of the formulas from Appendix B. Adding 
the results of all Yi integrations, we finally obtain the total contribution to the 
divergent part of / from all Yi terms (in Minkowski space): 
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r(4 - 2u) 



24 7^ 



2 2 

m — p 
n-p 



p 



p2 _ ^2 



p"^ +B 

p2 _ ^2 



+ 



N • ( 8, 22, -39, 89, ) + 



m / m 



p 



N • ( 0, -8, 4, -4, ) + 



m 



p. 



p^ 



- 1 In 1 - ^ - In(m^) 



N • ( 8, 0, 0, 0, 24 ) + 



2 2 

m — p 
B + p^ 



In 



— p"^ 



+ ln(S + m^) 



N - ( 0,-16, 16,-32, 24 ) 



, „ n-pn -p 

where B = — 



n 



2 ' 



L2iz) ^ 



^ Inxdx 
'o X- 1/z ' 



(5.9) 



the vector N is given by: N = ( m, ^, 



n-pi^* n*-pi^ p^ tji 
n-n* ' n-n* ' n-p 



(5.10) 



(b) Integrations over all parameter-space 

Finally, we must integrate the Yq terms over the finite parameters A, /3, G, 
b, h, and a (cf. Eqs. (5.1) and (4.2)). Excluding the r(2e) factor produced by the 
A integration, the integral of a Yq term will generally include a divergent term 
proportional to e~^, plus a finite term which is independent of e, plus terms 
which go to zero as e ^ 0. Due to the r(2e) factor, the divergent term will 
contribute to /_2 (cf. Eq. (4.8)), while the finite term will contribute to 
Since we need both /_2 and 7_i, we must evaluate both the divergent and finite 
parts of the integral of Yq. The divergent part comes from integration of Yq at 
its sub divergences, while the finite part comes from integration over the whole 
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region Since the subdivergences are part of the integral over the whole region, 
we can extract both the divergent part and the finite part by integrating Yq over 
aU of $. 

The functional form of the Yq terms is simpler than that of the Yi terms, 
because the former involve only the first term of the series in Eq. (4.9) and 
are, therefore, independent of H. However, when the denominator factor a is 
present in the original integrand f{q, k) (i.e., when G 7^ 0), integration of the Yq 
terms is still complicated by the presence of the factors = a(l — a) — G"^ and 
D_L = Xh — (see, for instance, Eq. (5.2a)). In the integration of the Yi terms, 
we simplified these factors by dropping all terms except those of lowest degree in 
{a,G} or {1 — a, G}. This procedure is, unfortunately, of no use here, because 
higher-order terms contribute to /_i via the finite part of the integral of Yq. Our 
new strategy is, therefore, as follows: we first simplify the and D± factors 
by changing the integration variables from {G, a. A, /3, b} to {V, U, X, r, t), 
according to the following definitions: 

G = Va, a=- — ^ , A = a(X + V2), /3 = ra , 6 = to. (5.11) 

Note that this change of variables does not involve the parameter /i, and 
is, therefore, inapplicable to any integral containing the factor n-k in the de- 
nominator of the original integrand f{q,k). This restriction is not as severe as 
it might seem, because we may always interchange q and k before integrating, 
as mentioned in Section 3(c). Hence the transformation (5.11) may be applied 
to all integrals except those with both n-k and n-q in the denominator of /. 
Fortunately, all such integrals of interest here (see Table 4) are harmless anyway, 
since they give rise to Yq terms which are proportional to n*n*, and hence vanish 
in the light-cone gauge. 
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Under the change of variables (5.11), the boundaries a = and a = 1 (at 
which all subdivergences occur) are transformed to C/ — > oo and U = V = 0, 
respectively, while the integrations in Eq. (4.1) are transformed according to: 

/•i i-i-G j-i- a j-X j-a 

dG da dX dp db, 

Jo Jg Jg Jg Jg 

i i ill (5.12) 

pi poo pU pX+V^ pi 

dV dU a'^ dX dr dt , 

Jo Jv-v^ Jv-v^ Jv Jv 

where d is the number of factors in the denominator of the original integrand 
/, a*^ is the Jacobian determinant of the transformation, and parameter integra- 
tions corresponding to absent denominator factors are to be omitted, as already 
explained in Section 4(a). 

Returning to our example, we now integrate Eq. (5.2a) over the new parame- 
ters. From Eqs. (5.11) and (3.7), we deduce D\\ — a?U and — a^X , and since 
J db does not appear in Eq. (4.2), J dt is absent from the transformed integral: 



$ 2D?,D\ ^ Jq Jv—v^ Jv—V^ Jv IU'^X^ ^ 



2 



dV / dU a^^ 
Jv-v^ 



+ 



1 + e (l + e)eC/2 



(5.13) 



Since the old parameter a is equal to the complicated function [U + + 1)"-*^, 
the new expressions a^U and a^X for D\\ and may not seem to represent 
a significant improvement over the old ones. However, when all factors of a in 
the transformed integrand, including the Jacobian a^, are collected, it turns out 
that the net power of a is equal to 2e in every case. This fact suggests that we 
may be able to obtain the divergent and finite parts of the Yq integrals by writing 
(t/ + + 1)~^^ in series form, and keeping only the first one or two terms. 
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In order to avoid the catastrophe of Eq. (2.11), we need series which converge 
uniformly over the whole domain of the U, V integrations. For U > 2, the binomial 
series 



26 TT-2e 



^ _ 2e {V^ + 1) 2e (2e + 1) {V^ + 1)^ _ 
1 t/ ^ 1 2 C/2 



(5.14a) 



is suitable, since its convergence accelerates with increasing U. On the other hand 
for C/ < 2, an exponential series is appropriate: 

{U + + 1)-^^ = 1 - 2e\ii{U + + 1) + 2e2(ln([/ + 1/2 + 1))^ - ... . 

(5.14b) 

In practice, one may avoid splitting up the integration region by using Eq. (5.14a) 
for integrals which diverge as t/ — > oo (e.g.: the first term in brackets in Eq. 
(5.13)), and Eq. (5.14b) for integrals which either diverge for finite U, or which 
do not diverge at all (e.g.: the other two terms). This strategy is valid, because 
in those parts of the integration region where an integral converges, the finite 
part of the integral may always be found by setting e = 0, either before or after 
integration. If we set e = in Eq. (5.14a) or (5.14b), all terms except for the 
first one go to zero, and the two series become identical. 

It turns out that, even at the points in the region of integration where a Yq 
integral diverges, only the first term of the appropriate series contributes to the 
divergent and finite parts of the integral. To see why, we recall from Section 4(b) 
that the degree of the integrand of every divergent finite-parameter integral is 
zero at the point in parameter-space where the divergence occurs; in other words, 
every integral diverges like f dx/x at a; = 0; if the degree were greater even 
by one, the integral would not diverge. In the series (5.14a), all terms but the 
first are of higher degree in 1/U than the factor (U + V'^ + 1)~^^, so that the 
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integrals involving these terms will remain finite as U — > oo. The accompanying 
factors of e in Eq. (5.14a) ensure that these terms do not even contribute to 
the finite part of the Yq integral as e ^ 0. Hence we may drop all terms of 
the series, except the first one. A similar argument applies to the series (5.14b), 
where ln{U + V'^ + 1) ^0 linearly as U,V ^ 0. 

The preceding analysis tells us that the complicated factor (U + V'^ + 1)~^^, 
appearing in the transformed integrands of the Yq integrals, may always be re- 
placed, either by U~^^ or by 1. Once this replacement has been made, integration 
of the Yq terms over t, r, X, and U is trivial, while final integration over V is 
easily accomplished with the help of the formula 

here C and D are linear functions of e. 

The results of the integration of all Yq terms are shown in Appendix C. When 
these results are multiplied by the appropriate coefficients and then added, the 
total contribution to the divergent part of / from all Yq terms reads as follows (in 
Minkowski space): 

^-•-^"■ C"W'°' ^t8,-4, n.,0)). (5.10) 

where N has already been defined in Eq. (5.10). To obtain the grand total of 
all divergent parts of /, we simply add Iyq from Eq. (5.16) to ly^ from Eq. (5.9). 
Thus (cf. Eq. (2.7)) 

= ig^T^DAT^DTicT%B[lYo + ly. + finite terms]. (5.17) 
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6. Concluding Remarks 

In this paper we have developed a powerful new integration technique for 
multi-loop integrals, called the matrix method, and used it to compute the di- 
vergent part of the overlapping two-loop fermion self-energy function zE in the 
light-cone gauge. The overlapping self-energy diagram (Fig. 1) is the most chal- 
lenging of the two-loop diagrams, especially in a noncovariant gauge such as the 
light-cone gauge. (The contributions from Figs. 2(a) to 2(d) are calculated in the 
sequel to this paper.) For completeness, we remind the reader of the one-loop 
result (Fig. 3) in the light-cone gauge [15]: 



The main results for the overlapping self-energy function (cf. Fig. 1) may be 
summarized thus: 

(i) The total divergent contribution, given by the sum of Eqs. (5.9) and (5.16), 
contains both simple and double poles. The double-pole term from Eq. (5.16) 
reads 



which is seen to be local, even off mass-shell. The contributions to [zE] 2 _ 2 



from results (6.1) and (6.2) are, therefore, strictly local. The coefficient of 
the single pole, on the other hand, has also non-local terms, as seen from 
Eq. (5.10). (Both types of poles are local on mass-shell.) The locality of 
the double-pole term (6.2) strongly suggests that the complete fermion- mass 
counterterm will likewise be local, but confirmation will depend on the results 
from Figures 2(a) and 2(b). 




-|- finite terms. (6.1) 




(6.2) 



'p =m' 
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(ii) The overlapping fermion self-energy integral contains a maximum of seven 
propagators. Detailed analysis of the singularity structure of the correspond- 
ing seven-parameter integrals in the "user-friendly" set S = {A; A, /3, G, b, 
h, a}, Eq. (3.4), reveals that the first simple pole originates from integration 
over the infinite parameter A, while the second simple pole (and double pole 
overall!) arises from integration over some of the finite parameters A, P, G, 
b, h, and a (subdivergences) . 

(iii) The matrix method is amazingly powerful, being applicable not only to mas- 
sive and massless integrals, but also to covariant integrals, as depicted in Ta- 
ble 3, and to noncovariant integrals, i.e., those containing the factors {k-n)~^ 
and/or {q-n)~^ (Tables 1, 2, 4, 5, and 6). The success of the technique de- 
rives, quite simply, from combining the 2a;-dimensional momentum vectors 
Qij, and /c^, and then integrating over 4c<;- dimensional Euclidean space in a 
single operation. It is this compact procedure which yields exact formulas 
(at an intermediate stage), whose analytic structure simplifies the ensuing 
parameter integrations tremendously. We note in passing that the matrix 
method works equally well for axial-type gauges, notably the temporal gauge 
(n^ > 0) and the pure axial gauge (n^ < 0). 

(iv) Although it has not been possible to corroborate our final result against 
other calculations (there just aren't any!), we have nevertheless been able to 
check all the covariant integrals listed in Table 11. And these agree exactly 
with the divergent parts of the double covariant integrals currently used in 
computing radiative corrections in the Standard Model [10,11]. Accordingly 
we feel reasonably confident that our final answer for the overlapping integral 
/ in Eq. (2.7) is correct. 
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Figure Captions 

Fig. 1 Two-loop overlapping quark self-energy function (solid lines denote quarks, 

wavy lines gluons). 
Fig. 2 Other two-loop diagrams for the quark self-energy. 
Fig. 3 One-loop quark self-energy function computed in ref. [15]. 
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Appendix A: Useful momentum-space integrals. 

The following formulas were derived from Eqs. (3.5) to (3.14). Note that 
z = {k^, q^, k3, q3, . . . 



J[k^] = r^J[l], 



(A.l) 



J[n*-qkJ = 



J[n* -qn-k ^] — 



J[n* -qn-k^] — 



J[n* -qq-k] 



where J[l] 



Gn* 



n -sr 



^ 2ADn 



n -s s 



2ADn 



n -sn-r^ + ) J[l], 



7r\4-2e Q-AH 
~A 



(e = 2-uj), 



(A.2) 
(A.3) 

(A.4) 



^ , /3n*-pi + 2Gnlif + f\\) \ ^^^^ 
n -sn-rf + — — ] AM, (A.5) 



(A.6) 



(A.7) 



(A.8) 



H= ^ f = (b + P-G){p'' + m^)-ibr + /3s)-p. 



(A.9) 



Since J is a linear functional, we can use J[q^] to deduce integrals such as J[n*-q] 
and J[^], etc. 
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Appendix B: Formulas for Yi integration 



Integration of Yi terms requires integration over finite parameters of ex- 
pressions involving Hq = Ha->o and Hi = h^i- -^o is given by Eqs. 
(5.5), which in turn were derived from Eq. (A. 9). We may similarly derive 
Hi = p| 6 (-R — 6 — bS) , where R and S are the same as in Eqs. (5.5). 

Formulas (B.l) to (B.ll) may be obtained by elementary means. Since = 
P\\ it follows from Eqs. (5.5) that p\{R - S -1) = rri^. 



dX 



InHn 



ln{m'^) - 2 + {R-S) In 



R-S 
R-S-1 



(B.l) 



db InHi = In(m^) -2 + 



R 



S+1 



In 



R 



R-S-1 



: / dbblnHi 
Jo 



In(m^) - 1 + 



R 



R-S-1 



R 



S+1 



(B.2) 



(B.3) 



1 rA /3pl 
dX / dp 

^0 



XHo 



S + 1 
R 



where L2{z) 



-U 



^ \n.x dx 
x-l/z 



(B.4) 



Formula (B.4) was derived with the help of the substitution (3 — yX , followed by 
integration over A and then integration by parts over y. The same substitution 
is helpful in the derivation of some of the integrals (B.5) to (B.ll) in Table 9. 



In formulas (B.5) to (B.ll), W = 



In(m^), 1, In 



R 



R-S- 1 



(R-S) In 



R-S 
R-S- 1 



S+1 
R 



S' 
R 



(B.12) 
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so that, for example, the right-hand side of formula (B.5) is 



In(m^) - 4 + {R-S)ln 



R-S 
R-S- 1 



+ RL2 



S+1 
R 



- RU 



'S_ 
R 



Table 9 



E 



1 rX 

dX df3 E 



^0 



In Hp 
X 



P In Hp 
A2 



(3 In Hp 
A 

X^Ho 

XHp 

2 

X^Hp 

XHp 



W 



w 



w 



w 



w 



w 



w 



1, 



1, 



1, S-ZR-2 



0, 



0, 



0, 



-1 



R 



- 1 



, 1 , i? 



R^ 



R 



-2^ 



, 3i? — S , 2i? 



' ^(5+1) ' S 



-R 



S+1 



R' 



s{s+iy 



1 , R 
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Appendix C: Results of Yq integration 



Each entry in the foUowing tables, when multiphed by 7r^'^r(4 — 2u!) times 
the Euchdean-vector expression at the right-hand side of its row of the table, 
gives the Yq portion of the divergent part of the integral of the expression at 
the left-hand side divided by the denominator at the bottom. Denominators are 
represented in the notation of definitions (3.15), and e = 2 — uj . 



Table 10 





-1 


-1 


+ i 
2€ 2 


+ i 
26 2 





n*/n-n* 

fi/ 




1 + 1 
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1 + 1 
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1-1 
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1-1 
t 


1-2 

t 


nl/n-n* 












i-1 






^ + 1 


i + 1 


1-1 


1 1 
26 2 


1-1 


iji* n-p/n-n* 




2 


6 


1-1 


1 + 1 









z 


D 




2 1 

e - 1 


U 


nk /fi 




1 


1,1 
26 ' 2 








7^n*-p / n-n* 


n-kJ/t 


11 


9 








-^||/8 




5 


3 








-A/8 




nq ^Kk 


n Q^Kk 


nqQ^K 
Table 11 


nqQA k 


nqQ Kk 




1 


1 + 1 


1 + 1 


1 + 1 


1 + 1 


2 
6 


1 


i 


1 + 1 


1 + 7 


1 + 2 


1 + 6 


8 
6 






q^Kk 


QAKk 


qQAK 
Table 12 


qQA k 


qQ Kk 




i 


lA 





2/6 


(4/6) -4 


7/in*-p/n-n* 




n AKk 


nq aK 


n Qa k 


n Q Kk 
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Appendix D: Computational summary of a typical overlapping integral. 



We summarize the main steps in the computation of the double integral, 
corresponding to the fifth row and second-last column of Table 1, Section 3(c): 



IeU] = f (f-q I (f-kf, with / 

J E J E 



nqQA k 



(D.l) 



Using definitions (3.15), we find that 



lE[f] = I I 
J E J E 



E n-q q"^ [{p — g)^ + m^] [{p — k — g)^ + m^] k"^ 



(D.2) 



1. First, we perform the exponential parametrization (2.9) with Schwinger pa- 
rameters CKi, ... , CKy (q!5 = Q!7 = 0); we then replace these parameters by 
the "user-friendly" set {A; X, l3,G,b, h, a}, defined in Eqs. (3.4). Thus (cf. 
Eq. (4.2)), 



Isif] = / dai ... j da^ / da^ J 
Jo Jo Jo 



-n*-q ^ 



Q!7 — > 




A^dA J 



^ Jo 



-n -q 



nt 



h ^ a 



(D.3) 



(D.4) 



with J defined by Eq. (3.3), and = dC da d\ j df3 . 

J^ Jo Jg Jg Jg 



2. Next, we perform the 4c<;-dimensional momentum integration in J with the 
help of formulas (A. 4) and (A. 8). In this way we obtain (cf. Eq. (4.4)): 



n 



4-26 




dA 



Jo 



^26 



-AH 



n*-s ^ + 



aiji* 



2ADn 



(D.5) 



with Dj^, r, s and H being defined in Eqs. (3.7), (3.14) and (A. 9). (In 
Eq. (D.5), and all subsequent equations, b ^ G, h ^ a, and e = 2 — a; .) 
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3. Integration over the infinite (type I) parameter A yields the first simple pole 
(cf. Eq. (4.6)): 

Hf] = / [r(2e) Jo + r(l + 2e) Ji ] , (D.6) 

with Jo and Ji given in Eqs. (4.7). The term with r(2e) diverges as e — > 
{uj —>■ 2). Additional subdivergences will arise from integration over $. 

4. Both Jo and Ji include factors of H~^^. According to Section 4(b), there 
are no subdivergences at the zeros of H; therefore, we use the exponential 
series for H~^^ given in Eq. (4.9) to get (cf. Eq. (5.1)): 

lE[f] = / [r(2e)ro + r(l + 2e)Yi +0(e)], (D.7) 

with Yq and Yi shown in Eqs. (5.2). As discussed in Section 5, we require 
both the finite and divergent parts of J^Yq, but only the divergent part of 

5(a). In this particular example, all subdivergences occur at a = G = 0. (In 
general, according to Section 4(b), subdivergences may also occur at a = 1, 
G = 0.) To find the divergent part of Yi, therefore, we integrate only the 
portion of Yi of least degree in {a, G} (cf. Eqs. (5.3) to (5.7)): 

Yi = dG da dX d(3[a^-'^E +0(a^-^)], (D.8) 
J* Jo Jg Jg Jg 

' dG f a^-^ da I dX I d(3 E + finite terms , (D.9) 

JG Jo Jo 

with E being defined in Eq. (5.4). Since E is independent of a and G, 

[ Yi = - [ dX [ dp Ee=o + finite terms. (D.IO) 
J$ ^ Jo Jo 
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5(b). Integration over A and (3 with the help of Eqs. (B.9), (B.8), and (B.5) yields: 



W 

e 



-1' 77^' 1, ^ ) - Tr(l: -4, 0, 1, R) 



S+1 



R 



S{S+1)' s 



+ finite terms , (D.ll) 



with W being defined in Eq. (B.12), and R and S in Eqs. (5.5). 

6(a). To integrate Yq over we first change the integration variables from {a, G, 
A, /3} to {U, V, X, r} via Eqs. (5.11). From Eqs. (3.7) we deduce D\\ = o?U 
and D_L = a^X, so that the integral of Eq. (5.2a) over $ becomes (cf. Eq. 
(5.13)): 



Yo = dV dU 
Jo Jv-v^ 



u 



a" I dX I dr 

v-v^ Jv 



* „-5+2e 



^* a 



dU a 



26 



V-V^ 



(y_y2)l+e y_y2 

TT7^ (l + e)et/2 ~ et/2-e 



, (D.12) 



with a= (U + V^ + 1)~\ according to Eqs. (5.11). 

6(b). As explained in Section 5(b), we obtain the divergent and finite parts of 
Yq by replacing a^^ with U~^^ in those terms whose integrals diverge as 
t/ — > oo, and with 1 in the other terms. Thus, 



Yq = i- I dV 



e(l + e) e(l + e) e(l - e) 



+ 0(e) , 



- 1 + 0(6) 



(D.13) 



The divergent part of IeU] may be obtained by combining Eqs. (D.7), (D.ll), 
and (D.13). 
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